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Abstract 

It lias recently been shown that the third law of thermodynamics is violated 
by an entire class of classical Hamiltonians in one dimension, over a finite 
volume of coupling-constant space, assuming only that certain elementary 
symmetries are exact, and that the interactions are finite-ranged. However, 
until now, only the existence of such Hamiltonians was known, while almost 
nothing was known of the nature of the couplings. Here we show how to 
define the subvolume of these Hamiltonians — a 'wedge' in a d-dimensional 
space — in terms of simple properties of a directed graph. We then give a 
simple expression for a specific Hamiltonian Ti* in this wedge, and show that 
7i* is a physically reasonable Hamiltonian, in the sense that its coupling 
constants lie within an envelope which decreases smoothly, as a function of 
the range I, to zero at I = r + 1, where r is the range of the interaction. 

I. INTRODUCTION 

It is sometimes statedS that all materials in their lowest-energy states are perfect crystals, 
i.e. that matter at zero temperature is characterized by periodic order of the atoms. If this 
is true, it follows that disorder persisting to low temperatures, in amorphous solids for 
example, must be interpreted as a result of trapping of the system in a metastable state. 
Another consequence is that the third law of thermodynamics holds in the Planck form,l'l 
which states that the entropy density tends to zero as the temperature T — 0. 

These statements, however plausible,! have not been proved. One might aim to prove 
that any physically reasonable microscopic Hamiltonian describing a material has a unique 
ground state which is spatially periodic. More generally, one wishes to know the minimal 
conditions on the Hamiltonian sufficient to guarantee a periodic ground state. Both problems 
are unsolved in general, but some progress has been made, particularly for one-dimensional 
systems.i Radin and Schulmani showed that if attention is restricted to a one-dimensional 
system of interacting classical units ('spins'), each of which can exist in a finite number k of 
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distinct states, with no interactions beyond a spatial range r, then there exists a ground state 
which is periodic with period at most k^. In particular, if the ground state is nondegenerate 
then it has perfect periodic order. 

Thus, in this class of model systems, disorder can occur only if the ground state is 
degenerate. This result can be strengthened by the observation that degenerate ground states 
occur only 'rarely', in the sense that they require fine-tuning of the system's parameters 
(coupling constants) to precise values.ili'Q In other words, degeneracy occurs only on a set 
of measure zero in the space of Hamiltonians. In the absence of accidental degeneracies, 
then, Radin and Schulman's result implies that the ground state of such a discrete classical 
system is always periodic. 

Recently, however, Canright and Watsoni (CW) have shown that this picture must be 
modified if the system is constrained by an exact symmetry. The idea that symmetry can 
imply degeneracy is familiar. For the discrete classical chain, CW showed that, under suit- 
able circumstances, the degeneracy arising from symmetry can result in a nonzero entropy 
density, throughout a finite volume of the space of coupling parameters. In this phase, 
termed a D-pair phase, almost all the ground state configurations are aperiodic. The D-pair 
phase is robust, in the sense that it is not sensitive to small perturbations in the coupling 
constants defining the Hamiltonian, as long as these perturbations respect the symmetry 
and the restriction to interactions of range r. It is also sufficiently robust to persist to finite 
temperatures. 

CW considered two symmetries in detail: spatial inversion (J), and spin inversion (S). 
They showed rigorously that, for S symmetry, D-pair phases exist if and only if k is odd, 
while for / symmetry, they exist for > 3 and r > 2. The Ising {k = 2) case is exceptional 
in that D-pair phases occur with I symmetry only for range r > 5. 

Although the CW proof is constructive, in the sense that it provides a method for finding 
all possible D-pair configurations for given k and r, there are immediate open questions. 
The CW result establishes existence or nonexistence of D-pairs in each case, but gives no 
information on the characteristics of the region in the phase diagram (the space of coupling 
parameters) occupied by the D-pair phase, when it exists — except that it has finite volume. 
One would like to know the size and location of the D-pair region. The location is important, 
since D-pair phases are of little interest unless they occur in a physically reasonable part of 
the phase diagram. For example, consider a Hamiltonian whose coupling constants increase 
with spatial separation, and then drop to zero beyond the cutoff range r. We consider such 
a Hamiltonian to be 'unphysical'. Conversely, if the D-pair region includes Hamiltonians 
whose couplings decrease smoothly as a function of interaction range /, reaching zero at 
/ = r + 1 (or before), then we would claim that the case has been made that physically 
reasonable Hamiltonians can give ground states violating the third law. 

In this paper we investigate these questions. After Sec. ||, which reviews the formalism 
used in the construction of D-pairs, we derive in Sec. ^ results which characterize the geom- 
etry of the D-pair region in terms of the combinatorial properties of the corresponding graph 
cycles. In Sec. |V|, we provide a simple construction for writing down an explicit Hamiltonian 
corresponding to any given D-pair. We prove that this Hamiltonian has couplings which 
fall off approximately linearly with distance, which shows that it is indeed possible to have 
D-pair phases without pathological Hamiltonians. 
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II. GRAPHS, CORRELATION POLYTOPES AND D-PAIRS 



The system of interest is composed of interacting classical units, forming an infinite one- 
dimensional chain. Each 'spin' cxj can take k distinct values, which we label 0,1, . . . , k — 1. 
A general Hamiltonian having interactions of maximum range r can be written 

-f^ = • • • (1) 

i 

where the sum is over all sites. Our interest is in ground states of H, which are those 
configurations {cTj} that minimize the energy density in the thermodynamic limit, i.e., 7i = 
H/N with the number of sites — > oo. In particular, we seek ground states which do not 
require fine-tuning of coupling parameters to precise values; hereafter we restrict the term 
ground state to mean minimum-energy states which are robust with respect to small changes 
in the Hamiltonian. (A more precise definition in this context is given in Ref. |^.) 

It is very useful to represent the Hamiltonian pictorially as a directed graph G^'^^ with 
energy weights assigned to the arcs.ili^ The graph has k^ nodes, each representing a possible 
sequence of r spins in the system. The arcs in the graph correspond to the operation 
of spatial translation in the chain by one unit: a directed arc connects two nodes if the 
rightmost r — 1 spins of one agree with the leftmost r — 1 spins of the other. The arc 
pointing from the node (ctq, . . . , cr,._i) to the node (ai, ct2, . . . , cXr) is assigned an energy 
weight f{crQ,ai, . . . ,0"^). Any spin configuration of the chain is represented by an infinite 
path in the graph, and its energy density equals the average weight of the arcs in the path: 
e = E/N = J2a /(c)^o-, where n„ is defined as the average occurrence of an arc cr in the path. 
Thus, each spin configuration is characterized by its arc densities {rio-}- The arc densities 
are not all independent, since they satisfy flow constraints! which state that at each node 
the sums of incoming and outgoing arc densities are equal. In addition, they satisfy the 
inequalities < < 1. 

In this language, the Radin-Schulman result is easily understood. Any path in a graph 
may be decomposed into simple cyclesi (SCs), where a SC is a closed path not visiting any 
node more than once. If G^l^^ has a unique SC with lower energy per spin than any other, 
then the nondegenerate periodic ground state of H is generated by repetition of that SC; 
if there are two or more lowest-weight SCs, then there is always a periodic ground state 
generated by repeating one of them. In either case, the period of the periodic ground state 
is at most the number of nodes in Gl!'\ which is k^. 

That these SCs are true ground states, in our restricted sense of being stable to perturba- 
tions in the Hamiltonian, is readily understood using the idea of the correlation polytopei'0 
P^{^\ The spin correlations are defined by 

^« = (^r<i...af;.), (2) 

where a denotes the sequence of integers {po,Pi, ■ ■ ■ ,Pr)', there are d = {k — l)k^ independent 
spin correlations, given by the values Pi = 0,1, k — 1 with pq ^ 0. To any configuration of 
the chain corresponds a dimensional vector s of correlations, and any Hamiltonian density 
can be written as a linear combination, Ti = —Yl, JaSa = —J ■ s, where the Jq, are the d 
independent coupling parameters. However, the mapping from configuration to correlation 
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vector is not one-to-one, and not all correlation vectors represent feasible configurations. 
Specifically, the correlations and the arc densities are linearly related (Sec. IV), and the 



constraints < rio- < 1 on arc densities translate to inequalities on the correlation vector. 
They constrain s to lie inside a convex polytope, and this is the correlation polytope P^^^^. 

Because the Hamiltonian is a linear function of the correlations, the ground states which 
are robust to small changes in couplings J a are precisely the vertices of P^'^\ By a simple 
argumentilB the vertices can be shown to be in one-to-one correspondence with the SCs of 
G'^\ and we arrive at the result that the ground states are 'almost always' periodic. One 
can enumerate all possible ground states by finding all SCs of the graph. 

The argument just sketched does not apply when a symmetry X is imposed, forcing 
symmetry-related arc weights to be equal. If the lowest-weight SC is not symmetry- invariant, 
there must be a pair of degenerate lowest-weight SCs. If these do not share a node, there 
exist two symmetry-broken periodic ground states. If they share one or more nodes, the 
domain wall energy between them is zero, so that there are infinitely many degenerate ground 
state configurations, most of which are mixtures with a nonzero density of domain walls. 
The latter case is the D-pair phase, so called since it comes from a pair of symmetry-broken 
configurations, and is characterized by Degeneracy (infinitely many ground states, yielding 
a nonzero entropy density) and Disorder (almost all the ground states have no long-range 
order) . 

A simple example0 illustrating the idea of a D-pair is the /c = 3, r = 1 model 

H = -{a^) + {alal,), (3) 

where the spins (jj can take the values and ±1 and the angular brackets denote an average 
over the chain. Ti is invariant under spin inversion (5) symmetry, a —a. It is useful 
to transform to the variables r = 2cr^ — 1, which take the values ±1. The Hamiltonian 
becomes, apart from irrelevant constants, Ti = (rjTj+i), the Ising antiferromagnet. Its 
antiferromagnetic ground state, when transformed back to a variables, is (. . . ±0±0 . . .), 
where each ± spin can take any value independent of all the others. 

The degeneracy and disorder in the ground state appears in this example as a trivial 
consequence of the double-valued transformation between r and a. What makes it special is 
the fact that these properties are stable to perturbations in the Hamiltonian, provided these 
respect the symmetry and the restriction to range 2 interactions. This follows from two 
facts, (i) The ground states are lower in energy density than any other (periodic) state by a 
discrete amount, so a sufficiently small perturbation cannot create a new ground state, (ii) 
The only allowed perturbation terms are those involving correlations already in 7i, plus the 
additional correlation (cTjCTj+i), and each of these takes the same value on all ground states, 
so the degeneracy is not split. This is the robustness characteristic of a D-pair. 



III. CHARACTERIZING THE D-PAIR REGION 



To study D-pairs for general k and r, CW introduced the concept of the reduced graph 
^G^^\ Not all symmetry- related pairs of SCs of G^^'> correspond to possible ground states in 
the presence of symmetry, because the equality of symmetry-related arc weights can imply 
the existence of a third SC with lower energy than the original pair. We refer to this situation 
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as decomposition of a SC pair. The definitions of the symmetry-reduced graph ^G^J^^ and 
its SCs are tailored to take care of decomposing SCs, in such a way that the possible ground 
states are in one-to-one correspondence with SCs of -^G^^ For S symmetry, the reduced 
graph is constructed by identifying each node or arc with its inverse, and SCs are defined 
as usual as paths which do not self-intersect. For / symmetry, the definition of the reduced 
graph and its SCs is more involved; we refer the reader to CW for the technical details, 
including the classification of SCs into four topological types. 

The reduced graph ^G^J^"* allows the enumeration of the ground state spin configurations 
for all D-pair phases with a given k and r. Here, we address the question of the region 
in the phase diagram in which a given D-pair phase is stable. By the phase diagram, we 
mean the d^^^^-dimensional space (reduced from d dimensions by the constraints arising from 
symmetry) of the coupling parameters J^. 

First, let us discuss the problem unconstrained by symmetry. We ask, what is the region, 
W , of J-space in which a given configuration (i.e. SC) u is the ground state? In principle, it is 
a region bounded by hyperplanes corresponding to the inequalities 'Hiyj) < Tiioj'), where u' 
ranges over all other SCs. However, in general some of these inequalities are redundant. We 
wish to determine the minimal set of inequalities needed to specify W fully. The following 
two lemmas provide a solution to this problem. 

Lemma 1. Suppose the SC uj corresponds to a vertex v of the correlation polytope. The 
region W of the phase diagram in which is a ground state is specified by the inequalities 
J- (v — v') > 0, where v' ranges over the vertices neighbouring v, i.e. those vertices connected 
m Pj'^) to V by a one- dimensional edge. Furthermore, this set of inequalities is minimal, in 
the sense that if any one of them is omitted the resulting region is strictly larger than W . 

Proof: Let {vi, V2, . . . , Vg} be the vertices of P^^"* with vi = v, and suppose {v2, V3, . . . , Vp} 
are the neighbours of vi. Then W is the set of J such that J • (vj — vi) < for 2 < i < q, 
and we define W to be the set of J such that J • (vj — vi) < for 2 < i < p. Since W C W, 
to prove W = W we must show W C W. 

It follows from the convexity of Pj'^-' that the set of vectors Vj — vi, 2 < i < p, from vi 
to its neighbours spans the full (i-dimensional space. In fact, for j > p, 

p 

1=2 

for some > 0, with at least two ctj 7^ 0.0 If J 6 W, taking its dot product with both 
sides of (ID yields J ■ (v^ — vi) < 0, and hence J G W, as required. 

Define W" as for W but omitting one neighbour, say V2. Since V2 is a neighbour of vi 
and since Pj*^-* is convex, there exists a hyperplane of dimension d — 1 intersecting Pj'^) only 
in the edge joining vi and V2. Let J be a perpendicular vector to this plane from the origin. 
The sign of J can be chosen so that J ■ (vj — Vi) < for i > 2, and thus J G W", while 
J ■ (v2 - vi) = implies J ^W. ^ 

Lemma 2. Two vertices in Pj'^^ are neighbours if and only if the corresponding SCs of G^*^^ 
have zero or one contacts, where a contact is a consecutive sequence of one or more shared 
nodes. 
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(k) 

FIG. 1. Schematic picture of simple cycles in Gf- having two contacts. The two circles repre- 
sent the cycles, made up of arc sequences labelled with weights ni to n^, and the large dots are 
the contacts, which may consist of more than one node. 

Proof: The vector of arc densities, n, has dimension equal to the number of arcs, but the flow 
constraints (Sec. 0) constrain it to lie in a d-dimensional subspace which we denote P'. It is 
the image of Pj'^-' under a nonsingular linear transformation M from s to n (see Sec. IV for 



explicit relations). It follows that neighbouring vertices of P^^^-* correspond to neighbouring 
vertices of P'. Two vertices Vi and V2 are neighbours if and only if any point AiVi + A2V2 
(with Ai + A2 = 1) on the line segment joining them cannot be written as a weighted average 
of vertices in any other way. Suppose two SCs have two contacts, as illustrated schematically 
m Fig. I. Reco gnizing that the four arc sequences define four distinct SCs, we can consider 
a general convex combination including coefficients A3 for the inner cycle (arcs 2 and 3) and 
A4 for the outer cycle (arcs 1 and 4). A point on the line segment joining the corresponding 
vertices in P' has densities rii = = X and n2 = = 1 — X. Clearly, there are many 
convex combinations of vertices yielding the same densities; for example, if A < 1/2 we can 
take Ai = 0, A2 = 1 — 2A and A3 = A4 = A. Hence the two SCs correspond to vertices which 
are not neighbours. Conversely, if the SCs have fewer than two contacts there is only one 
way to express points in density space on the segment joining them as a convex combination 
of SCs, and so they correspond to neighbouring vertices. | 

We remark that the reasoning in Lemma 2 is very similar to that leading to decomposition 
of pairs of SCs (except that in Lemma 2 the SCs need not be related by symmetry). This has 
a simple geometrical interpretation. A symmetry imposes linear constraints on the coupling 
parameters J^, which means that the relevant space of correlations is a reduced polytope 
^P^'^^ obtained by symmetry projection of Pj^^\ When symmetry- related pairs of vertices of 
Pjf^^ are projected, the ones which become vertices of ^P^^^ are those which are connected 
by an edge. Indeed, the CW definition of SCs of ^G^*^^ (for X = S* or J) is constructed so as 
to include all cycles in G^*^-* which have at most one contact with their symmetry-partners. 
There is, however, one category of non-decomposing SC pairs which does not correspond 
to neighbouring vertices. It has the form of Fig. |l|, in the case that symmetry forces the 
weights to satisfy Wi = W2 and = 104. This situation occurs with / symmetry for a type 
four SC, when the two contacts are symmetry-inverses of each other. Because of the weight 
constraints all four cycles in the diagram have equal energy, and the original pair does not 
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decompose. It corresponds to two non- neighbouring vertices of a quadrilateral face of P^''\ 
such that all four vertices of the face map to the same vertex of ^Pj'^) under the symmetry. 

Lemmas 1 and 2 show how to determine the region in coupling space corresponding to 
a given ground state configuration; in fact, they provide an algorithm for doing this. From 
Lemma 1, only neighbouring vertices need be considered. Lemma 2 translates the concept 
of neighbouring vertices into properties of graph cycles. In terms of spin configurations, a 
contact between SCs means a common string of r or more spins. 

Let us now consider the analogous problem in the presence of a symmetry X. Since 
Lemma 1 relies only on the Hamiltonian density being a scalar product, it applies directly 
to the symmetry-constrained problem, i.e. the inequalities defining the stable region W come 
from neighbouring vertices of the reduced (projected) correlation polytope ^P^^\ Lemma 2 
also goes through unchanged in the case of S symmetry, since the definition of SCs for ^G[^^ 
is the same as that for G^J^\ However, Lemma 2 does not apply when the symmetry is /. 

As in the proof of Lemma 2, it is clear that a pair of SCs in ^G^^'^ represent neighbouring 
vertices in '^Pj*^-' if and only if there do not exist two or more new SCs of ^G'^'^ using only arcs 
from the original pair. Because SCs for I symmetry may, when unfolded into G^^\ represent 
pairs of intersecting cycles, there is more freedom to form these new SCs than in the absence 
of symmetry. We find that when the intersection does not contain a symmetric node, then 
one contact between the original SCs may be enough to imply new SCs. Specifically, we find 
the following: 

Lemma 2'. For / symmetry, two SCs of ^G^^^ correspond to neighbouring vertices if and 
only if one of the following conditions is satisfied: (i) they have no contacts; (ii) they have 
one contact and one of them is of type one (i.e. unfolds to nonintersecting cycles); (iii) they 
have one contact which includes a symmetric node. (In the last case, both SCs must be type 
two or three.) 

Let us illustrate our results with the example of A; = 3 and r = 2 for both S and / 
symmetries. As in Sec. ^ we shall take the allowed spin values to be cr = and ±1. 

Fig. H shows the graph '^G2^'*. It has 14 arcs, and 5 nodes each of which implies a flow 
constraint, leaving 9 independent arc densities, i.e. S-^'^ = 9. The 9 symmetry-invariant 
correlations are si = (cx^), S2 = (aiCTi+i), S3 = {aiai+2), S4 = (crfaf+i), S5 = (a.^af+a), plus 
four correlations involving three spins; the Hamiltonian density is written in terms of its 9 
coupling parameters as H = —J2 JaSa- The graph has 19 distinct SCs, 5 of which are D- 
pairs using the invariant node (00). For example, let us consider the D-pair SC uj = (00±). 
To find its stable region, we need only consider the 10 SCs which represent neighbours of 
u, according to Lemma 2. One neighbouring SC is the ferromagnetic state (00); comparing 
its energy to that of lj yields the condition Ji > 0, so we may set Ji = 1. For each of the 9 
other neighbouring SCs one can write down the corresponding inequality on the 8 remaining 
couplings directly from the graph. We do not list them here; let us merely display a typical 
solution: 

^ = -(^') + (^Mi) + KV+2)- (5) 

As in the example of Sec. || it is informative to perform the transformation r = 2cr^ — 1 to 
Ising spins. The Hamiltonian becomes 
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FIG. 2. The reduced graph ^G^^^ 



H = (TiTi+i) + {TiTi+2) + 2(r), 



(6) 



which represents an Ising model with antiferromagnetic nearest and next-nearest neighbour 
interactions and a magnetic field favouring the — state. It is easy to check that the ground 

state is ( h), or in a variables, (00±), as required. Thus we have the degeneracy and 

disorder characteristic of a D-pair. Its robustness to perturbations in the Hamiltonian follows 
from the fact that all 9 correlations take identical values on every degenerate configuration; 
i.e. no perturbation made up of si to sg can spht the degeneracy. 

The corresponding reduced graph for I symmetry is shown in Fig. ^. For ease in drawing 
we have distorted the symmetry line X into a circle and omitted the ferromagnetic arcs 
joining each /-invariant node to X. This graph has 32 distinct SCs, most of which are of 
type three, unfolding to invariant cycles in G2 . There are three (type two) D-pairs, of 
which we shall consider the example (00 — |-)/(00H — ). It has 15 neighbours according to 
Lemma 2', namely the ferromagnetic SCs, the type one SC (OH — )/(0 — h), and the 10 type 
three SCs which use the symmetric node (00). Thus there are 15 inequalities constraining 
the 14 distinct /-invariant correlations. Again, we do not list them here, but simply display 
a particular solution, which happens to involve only S'-invariant pairwise interactions: 



The ground state is (. . . 00±=FOO±=pOO . . .), where the spins in each (±=1=) segment may be 
chosen independently to be (H — ) or (— +)• Once again, one can check that this degeneracy 
is not split by any of the 14 possible /-invariant perturbing terms that may be added to the 
Hamiltonian. 



(7) 
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FIG. 3. The reduced graph ^G)^' . 

IV. EXPLICIT D-PAIR HAMILTONIAN 

The techniques described in the previous section can be used to find the region of stabihty 
corresponding to any given D-pair, for S or I symmetry. However, the analysis becomes 
tedious for large k and r. In this section, we describe a simple construction for finding a 
single point, H*, in the stable region. 

The constuction is based on the observation that the arc weights, which determine the 
Hamiltonian, can all be chosen independently, i.e. for any choice of arc weights there exists 
a corresponding Hamiltonian. (This statement should not be confused with the fact that 
the arc densities are not independent because of the fiow constraints, and hence that the 
map from arc weights to couplings Jq, is not one-to-one.) Given a D-pair defined by a SC u 
of , we define Ti.* as the Hamiltonian corresponding to the following assignment of arc 
weights: Wr is —1 if the arc r occurs in the unfolding of the D-pair SC into two intersecting 
cycles in Gl!'\ and otherwise. In terms of arc densities, Ti.* is 

n* = -Y.{n, + n,), (8) 

where the overbar denotes symmetry inversion. By construction, the given D-pair phase is 
the ground state of Ti.*: since the energy is the average arc weight, any of the D-pair spin 
configurations has energy —1, while other configurations use some weight arcs and have 
higher energy. 

In (H), Ti* is written in terms of arc densities. The latter are related to the correlations 
as follows. If T = {tq,Ti, . . . ,Tr), then 

The Kronecker delta, as a function of a spin variable a, can be written as a degree k — 1 
polynomial according to 
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fc— 1 

5<.. = - - = E Xpr,^'. (10) 

The product is over all spin values rj' not equal to 77, and the second equality defines the 
numbers Xpn as the coefficients in the polynomial expansion of the product. When this is 
substituted into the arc densities are given as a linear combination of spin correlations. 
From (D), this yields an expression for H.* in terms of the correlations. As an example, when 
this procedure is applied to the D-pair (±0) for = 3, r = 1 with S symmetry, the result is 
the Hamiltonian (H), discussed in Sec. |I|. 

Let us investigate further the structure of Ti*. Its expression in terms of correlations, 



from substituting (|T^) and into (||), is 



p k—l fc— 1 fc— 1 

^* = - E E E • • • E Xpor,Xpir,+i • • • Xp.r,+,(cTfVfji . . . df;,) + S.l. (ll) 

i = lPO = Opi=0 pr=0 

Here, Tj denotes the jth spin (using any arbitrary starting point) of the configuration defined 
by the SC cj, and p is the period of uj. The second term, not written explicitly, is the 
symmetry inverse of the first — that is, every r 1— > r. Consider a range / correlation s„ = 
{al'^alXi ■ ■ ■ (^f+i)^ where go and qi are both nonzero. If / < r, there are (r — / + 1) terms in 
(|n|) contributing to Jq,, the coupling parameter multiplying Sa in Ti*. We find 

p r—l 

J a — E/ E/ XOTj ■ ■ ■ XOTj + rn-lXqoTj + m ■ ■ ■ XqiTj + rn+lXOTj + rn + l + 1 ' ' ' XOrj + r ~^ i^^) 

j=l m=0 

The structure of this expression is seen most clearly if we consider initially the case of Ising 
spins, k = 2. Taking the allowed spin values to be a = ±1, we have Xo^ = 1/2 and Xir? = ^7/2, 
and we arrive at the result 

Ja = 2-(^+^) (r - / + 1) £ r/li . . . r][, + s.i. (13) 

i=i 

This result is quite significant. It says that the value of Ja is, apart from a constant 2~^'^~^^^ 
that we shall ignore, equal to (r — / + l)[ta + s.i.], where ta is the correlation Sa evaluated 
in the spin configuration {rj} of the D-pair. 

Since a correlation for Ising spins has magnitude at most 1, this implies the bound 
\Ja\ < (r — Z + 1). Further, we note that is expected not to be strongly dependent on 
the range / of the coupling Jq,. Roughly speaking, spin configurations that have long-range 
correlations tend also to be correlated at short range. This idea is borne out by explicit 

computations; for instance, for the /-symmetry r = 5 Ising D-pair (H hH — ), there 

are 23 symmetric correlations, each of which takes one of the values —1/7 or 3/7, with no 
systematic dependence on I. Thus, the dominant contribution to Ja comes from the factor 
[r — I + 1). As a function of distance, this represents a linear decrease to zero at the cutoff 
range I = r + 1. 

For k > 2, the situation is similar. Of course, the values of the couplings depend on the 
choice of the set of allowed spin values, which has been left arbitrary so far. However, for any 
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k there exists a choice with the property that xov is independent of 77, namely, letting {ex} 
be the complex fcth roots of unity. For this choice, the sum (0) simplifies as it did for Ising 
spins, yielding (apart from a constant) Ja = {r — I + l)[tQ, + s.i.]. Under the assumption that 
the correlations depend weakly on we find again that the dominant distance dependence 
of Ja is a linear fall-off to zero beyond the cutoff range. 



V. DISCUSSION 

The existence of D-pair phases is interesting from a theoretical point of view. However, 
we are not aware of any obvious candidate material for their realization in nature. We note 
that there is an entire class of materials, namely, layered solids or polytypes, which are well 
modelled by effective Hamiltonians such as those studied here. This class of materials is 
however quite large; and the few effective Hamiltonians that are known from this class do 
not show promise of having a D-pair phase as the ground state. (See the discussion and 
references in Ref. |lT].) 

Thus, there are significant obstacles to finding D-pairs in practice. However, the re- 



sults of Sec. 1^ of this paper remove one potential obstacle: the possibility that the only 
Hamiltonians exhibiting D-pair phases are pathological in the dependence of their coupling 
parameters on distance. We would like the couplings to decrease smoothly to zero at the 
cutoff range, otherwise it would seem unphysical to impose a rigid cutoff beyond which 
there are no interactions. We have constructed an explicit Hamiltonian for arbitrary k and 
r, which has D-pair ground states. Encouragingly, its couplings are very well behaved: as a 
function of distance, they fall linearly to zero at the cutoff. 

Remaining obstacles concern the robustness of D-pair behaviour. Although D-pairs 
are not destroyed by symmetric perturbations of sufficiently short range or by nonzero 
temperature, they are in general destroyed by including interactions beyond the cutoff. They 
are also destroyed by deviations from perfect symmetry caused, for example, by external 
fields — which may or may not be strictly zero, depending on the symmetry in question, 
and on the physical identity of the 'spins'. However, even when the degeneracy is broken 
in such ways, behaviour characteristic of D-pairs may be observable at a suitable energy 
scale. If the perturbations breaking the D-pair symmetry are small, they are not manifest 
except at very small temperatures. At low but nonzero temperatures, one would still expect 
to observe disordered states and nonzero entropy density. In that case, the techniques of 
Sec. of this paper apply directly to the problem of characterizing the D-pair region. 
(Possible experimental signatures of D-pairs have been investigated by Yi and Canright. 0) 

Another potential obstacle is the limitation to problems involving classical, discrete units. 
However, such models are likely to be good approximations for certain problems, such as 
stacking polytypes of crystals (see CW and references therein) where the 'spin' represents 
the discrete set of possible configurations of a single lattice plane. Another limitation is the 
restriction to one-dimensional models. The question of whether similar behaviour is possible 
in higher dimensions is unexplored, although certain frustrated two-dimensional models are 
known to have degenerate ground states of large periodicity.lll Finally, although we have 
shown that D-pair phases are possible with Hamiltonians that are not obviously unphys- 
ical, there may be more subtle physical reasons — arising, say, from quantum-mechanical 
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considerations — which may argue against effective classical Hamiltonians having D-pairs as 
ground states. For example, effective classical Hamiltonians representing the binding energy 
of mobile electrons in an ionic background tend to favour periodic ionic arrangements. H We 
leave these questions for future work. 
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